Non-equilibrium theory of charge qubit decoherence in the quantum point contact 

measurement 



Ming-Tsung Lee° and Wei-Min Zhang"^^ 

'^National Center for Theoretical Science, Tainan, Taiwan 70101, R.O.C. 
^Department of Physics and Center for Quantum Information Science, 
National Cheng Kung University, Tainan, Taiwan 70101, R.O.C. 
(Dated: April 19, 2008) 

A non-equilibrium theory describing the charge qubit dynamics measured by a quantum point 
contact is developed based on Schwinger-Keldysh's approach. Using the real-time diagram tech- 
nique, we derive the master equation to all orders in perturbation expansions. The non-Markovian 
processes in the qubit dynamics is naturally taken into account. The qubit decoherence, in partic- 
ular, the influence of the tunneling-electron fluctuation in the quantum point contact with a longer 
time correlation, is studied in the framework. We consider the Lorentzian-type spectral density to 
characterize the channel mixture of the electron tunneling processes induced by the measurement 
and determine the correlation time scale of the tunneling-electron fluctuation. The result shows that 
as the quantum point contact is casted with a narrower profile of the spectral density, tunneling 
electrons can propagate with a longer time correlation and lead to the non-Markovian processes of 
the qubit dynamics. The qubit electron in the charge qubit will be driven coherently. The quantum 
point contact measurement with the minimum deviation of the electron tunneling processes prevents 
the qubit state from the decoherence. 

PACS numbers: 03.65.Yz,73.23.-b,05.70.Ln,03.67.Lx 



I. INTRODUCTION 

With the rapid progress in nano-technology, the in- 
vestigation of quantum processes in systems couphng to 
a mesoscopical measurement device becomes a very ac- 
tive research field in recent years. P, [3, H, 0] Not only 
its practical application to quantum communication and 
quantum computation, 5] but also the theoretical in- 
terests in the measurement-induced quantum decoher- 
ence [3, @ have attracted much attention. In particu- 
lar, in the investigation of the solid-state quantum com- 
puter with charge qubits, the quantum point contact 
(QPC) has been served as an ultrasensitive electrome- 
ter [E 0j H, H, ■ the literature, the measurement of 
charge qubit through the QPC has been treated based on 
the so-called Markovian approximation. (Til. [T^ ITsl. IT3| in 
which the time scale of the qubit dynamics is assumed 
much larger than that of the tunneling-electron correla- 
tion in the QPC. This theoretical treatment based on the 
Markovian approximation describes only the qubit dy- 
namics in the time-asymptotic quasi-equilibrium state. 
Mesoscopically, the qubit decoherence occurs in the time 
scale of the same order of the tunneling-electron correla- 
tion time in the QPC, where the non-Markovian dynam- 
ics of the qubit is significant. Thus, a non-equilibrium 
description to the qubit dynamics is more desired. Re- 
cently, the considerations of the solid-state system pro- 
cessing in the non-Markovian regime have indeed received 
more attention. For instance, a local electron spin cou- 
pling with a nuclear spin bath through the Fermi contact 
hyperfinc interaction in which the electron spin dynam- 
ics is in the time scale shorter than the nuclear dipolc- 
dipole correlation time[l5|, and electron transports of in- 
teracting electron systems most likely involve the non- 



Markovian dynamics [T^ [l7| . In our previous workfT^. 
the non-equilibrium effect of the QPC on the qubit dy- 
namics was studied by treating the fermi energy fluctua- 
tion of the QPC reservoirs perturbatively. In this paper, 
a fully non-equilibrium theory describing the qubit de- 
coherence by the QPC measurement is developed using 
Schwinger-Keldysh's approach (isj. 

Historically, the Schwinger-Keldysh's approach is well 
developed to systematically treat the non-equilibrium 
dynamics of a many-body system. [l^ [2^ [2l[ For elec- 
tron transports in nano-devices, this approach has been 
used to study the transport current fluctuation and the 
full counting statistics in the single-electron transisto r 12 1[ 
l22j . the current fluctuation in Kondo system [23l. [2I. |25||. 
and the noise spectrum in the QPC measurement of the 
charge qubit [2f|. To treat the non-equilibrium effect 
of the electrical reservoir, literaturely the real-time di- 
agrammatic technique was constructed to diagrammati- 
cally calculate correlation functions of the electrical reser- 
voir order by order. [U [23| For the investigation of the 
non-Markovian dynamics in our system, an alternative 
real-time diagrammatic technique is developed. The 
master equation for the charge qubit dynamics is derived 
and expressed in terms of all orders irreducible diagrams 
to all orders. The non-Markovian processes in the qubit 
dynamics can be fully taken into account. The effect of 
the time variated reservoir fluctuation on the qubit dy- 
namics can then be explicitly studied in this formulism. 

In addition, the assumption of the constant hopping 
amplitude of tunneling electrons across the QPC barrier 
together with a constant density of state of the QPC 
reservoirs is usually utilized to specify the QPC struc- 
ture related to the two reservoirs band structure cou- 
pling to the two metal gates. It eventually leads to the 



2 



qubit dynamics in the Markovian hmit.[Tll. [T^ . [T3 | 
However, the QPC structure determines the correlation 
time scale of the tunneling electron fluctuation in the 
QPC. The non-Markovian processes of the qubit dynam- 
ics can emerge if a particular design of the QPC struc- 
ture is taken into account. We will consider in this work 
a Lorentzian-type spectral density to characterize the 
energy-level dependence of the hopping amplitude and 
the density of state. A close relation between the qubit 
decoherence and the time correlation of the tunneling- 
electron fluctuation shows that the qubit decoherence can 
be controlled through the measurement operation itself. 

We organize the paper as follows: The theory of the 
charge qubit measurement is presented in Sec. II. The 
real-time diagrammatic technique based on Schwinger- 
Keldysh's approach are developed in Sec. Ill, where we 
also derive the master equation for the reduced density 
operator of the charge qubit. In Sec IV, the qubit dy- 
namics is studied based on the master equation. The 
influence of random electron-tunneling processes on the 
qubit decoherence is illustrated in this section. Finally, 
a summary is given in Sec V. 



II. CHARGE QUBIT MEASUREMENT 

The charge qubit measurement using QPC is studied 
in the tunnel junction regimefl [13, [HI, [SM H HI . In 

this regime, the transmissions of all tunneling channels 
cross the QPC barrier are small enough such that the 
electron tunneling becomes sensitive to the qubit state. 
The qubit information can then be extracted from the 
output signal of the QPC, while the backreaction of the 
measurement to the qubit states is expected to be min- 
imum. The Hamiltonian of the whole system is given 

byiEIBEEHIiii 



H = Hs + Hb + H', 



(1) 



where Hs denotes the Hamiltonian of the charge qubit, 
Hb the Hamiltonian of the QPC with the electrical reser- 
voirs consisting of the source indexed by I and the drain 
by r. The charge qubit state is measured through the 
electron tunneling across the source and the drain. H' is 
the interaction Hamiltonian describing the electron tun- 
neling processes through the QPC with a qubit-state de- 
pendent hopping amplitude qri- Explicitly, 



Hs 



Hb 



eiai ai 



H' = '^{qrio^ai + qirolar). 



(2) 



(3) 



(4) 



rl 



We shall formulate the non-equilibrium theory of the 
electron tunnelings coupled with the qubit dynamics in 



the electron-hole representation, in which Hb and H' can 
be written equivalently as 



Hb = Y.{elatai+e1f3+(ii) 

l>0 



(5) 



r>0 



H' = ^ qri{a-+ai + /3r/3+ + /S^az 



+a+/3+) -f H.c, 



(6) 



where af^{ai^r) and /3/'^(A,r) are respectively the cre- 
ation (annihilation) operators of the electron and hole, 
the corresponding electron and hole energies ^ = e;,r — 

IJ-L,R and e\\ — ^l,r — ^i,r with respect to the chemical 
potential /iL of the source and /ir of the drain. Since the 
coupling Qri depends on the qubit state, it is indeed a 
coupling function of the qubit operator. Meanwhile, qn 
also depends on the measurement device structure. We 
remain the discussion of the QPC measurement to the 
qubit decoherence with a practical qri later. 

The qubit dynamics is determined by the master equa- 
tion for the reduced density operator psit) ^Tr^ [ptot{t)]: 
where ptot{t) is the total density operator of the whole 
system, and the partial trace Tr^ integrates over all the 
degrees of freedom of the QPC reservoir. From the Liou- 
ville equation ^Ptot = —i[H^ptot] for the total density 
operator, it can be shown that the reduced density op- 
erator ps{t) obeys the following equation of motion (the 
master equation) [13, 113, HI 

^os{t) = -^ [Hs,Ps{t)] - iTrs U'QH„(i - h) 
at L 



X Ptotih) 



dTK{t-T)*ps{T), (7) 



where K{t ~ r) * ps{t) =11x b\L' Q H,(t - t)L' pt,t(T% 
L' EE [H',], QHoit) = e-**[^s+ffB,]^ and to is the initial 
time that the interaction between the qubit and the QPC 
measurement turns on. The derivation of Eq. ([7]) can 
be found in Appendix A. The second term in the above 
equation has no contribution due to the particle number 
conservation in electron tunneling processes. Since the 
QPC output signal records the information of the qubit 
state through the interaction between the QPC and the 
qubit, the induced backreaction from the fluctuant reser- 
voirs will result in qubit decoherence which is in princi- 
ple non-Markovian. The non-Markovian dynamics of the 
qubit is described by the term — J^*^ dTK{t — r) * ps{t), 
which contains all influences of the fluctuant reservoir to 
the qubit dynamics. The qubit decoherence is thus to- 
tally governed by Eq. ([7]). In the following discussion, 
we will derive the master equation for the reduced den- 
sity operator in terms of a diagrammatic perturbation 
expansion of Eq. ([7]) . 
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III. REAL-TIME DIAGRAMMATIC 
TECHNIQUE 

We begin with the Schwinger-Keldysh's approach [l8l| 
to explicitly calculate the term — J^^ dTK{t — r) * ps(t) 
in Eq. ([7]). In the interaction picture, the interaction 
Hamiltonian Eq. ([6]) can be rewritten as 

Hi = J2 {<irimnt)Mt) + q+im+{t)Mt)}, (s) 

rl>0 

where the field operator ipk{t) is defined by 

V3fc(0 = afce-"^ + (9) 

with the index k — l,r labeled respectively for the source 
and the drain. In addition, K{t — t) * Ps{t) in the inter- 
action picture can be rewritten as 



K{t-T)*ps{T) ^Qs{t-to) 



-Trs 
-Trs 



Tp [e-'^l^o ''"'^ptotito)} H'^H'^ ) , (10) 



where Qs{t) = e j^j^j j^^^ denotes the closed- 

time-path integral along the closed-time-path contour p 
[isl . [igl [20l . [2l| with the range of the real time axis from 
to to t. The closed time path contains the positive branch 
and the negative branch. The positive branch coincides 
with the real time axis, and the negative branch is re- 
versed with respect to the real time axis. Taking pertur- 
bation expansion of Eq. JIT 



K{t -t)* ps{t) = K^°\t - t) * ps{t) 

+ K^^\t - t) * ps{t) + ■ ■ ■ (11) 

through the expansion [l9| 

Tp je""-^* ''''^=ptot(io)| = X! — T~ / '^^i ■ ■ ■ / 

n=0 ''P •'P 

xTp[H',^---Hi^Ptot{to)}, (12) 
we can diagrammatically illustrate Eq. (llOp as 



with 





















1 i-L-i^ 




fi't A\ 



hIh7\ 



Ptot (r) = Tp [e-'il^o '^^^'-p,,, Ito)] . (13) 



According to the Wick's theorem, all higher order cor- 
relation functions of the fermion field operators can be 
built up in terms of the unperturbed Green's function, [l^, 
[20l . [2l| The reservoir contour-ordered Green's function is 
defined by 

Gfc(i,t') = -*(rp(v'fc(t)^+(0)), (14) 



where the contour-ordering operator Tp orders the oper- 



ator ipk{t)i^'^{t') according to their time arguments along 



the p, and (.) denotes the non-equilibrium statistical av- 
erage. Explicitly, Gk(t,t') contains four components: 



(15) 



i.e. the time-ordered Green's function GF,k, the anti- 
time-ordered Green's function Gp^, the correlation 
Green's functions G<,fc(t,t') and G>_fc, defined respec- 
tively as 



(16) 
(17) 
(18) 
(19) 



GF.u{t,t') = -il^TiMmtit')) 
Gp^^it,t') = -i(T{Mmtit')) 
G^Mt,t')=i(i;+{t')Mt)), 
G>,fc(t,0 = -»(^feW^fc+(t')), 



where T is the normal time-ordering operator and T the 
anti-time-ordering operator. The unperturbed reservoir 
contour-ordered Green's functions are easy calculated. 



iGf^it,t') = Trs ^P+it')Mt)p 



= (l-./(4))e-^^'~*'' 

+/(e^)e-S(*-*'), (20) 



= (l-/(^D)e-^^^^*-''^ 

+f^,hy4it~t')^ (21) 

G^°i(i,t') =^(i-i')<i(i,t') + ^(i'-i)<l(i,t'), (22) 

G(,°^(i,t') =0(i-i')<U^,O+e(i'-i)Gl,°yt,t'), (23) 

Here, the electrical reservoirs are assumed to be in the 
thermal equilibrium state initially (at to), and /(e) 
is the Fermi-Dirac distribution function. These Green's 
functions describe the contractions of the unperturbed 
fermion pair with four different time contour-orderings. 



?,+ ^+'^ JO) 
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A. The diagrammatic rules 

In order to systematically trace out the reservoir de- 
grees of freedom to all orders in the perturbation expan- 
sion with a correct operator ordering for the remaining 
degrees of freedom of the qubit, we shall use the real-time 
diagrammatic expansion with the diagrammatic rules de- 
fined as follows: According to the interaction Hamilto- 
nian of Eq. ([8]), two kinds of the tunnelings, the forward 
and the backward tunneling, across the QPC barrier are 
involved. The diagrammatic representation of the for- 
ward tunneling vertex qriip^ipi and the backward tun- 
neling vertex are depicted in Fig. [1] (a). For 
the forward tunneling vertex, the incoming dashed line 
labeled by I represents the electrons (or holes) above (be- 
low) the chemical potential /i^, are destroyed (created) in 
the source, and the outgoing solid line labeled by r repre- 
sents the electrons (or holes) with the energy level above 
(below) the chemical potential are created (destroyed) 
in the drain. The QPC current due to this tunneling ef- 
fectively flows from the source to the drain. The coupling 
operator q^i is presented in the vertex with a filled circle 
for the forward tunneling. Similarly, the interaction as- 
sociated with the backward tunneling is depicted by the 
incoming solid line labeled by r and the outgoing dashed 
line labeled by Z, and the coupling operator presented 
in the vertex is denoted with a hollow circle. 



(a) 



forward tunneling backward tunneling 

ST TZ -. 1-^ 



^ 



I ~ 
I P 



is) 



(b) 



^y,G|') (.,,.,) 




FIG. 1: The diagrammatic representation, (a) The forward 
tunneling vertex qri{t)''i'r {'t)''Pi{'t) and the backward tunneling 
vertex qir{t)'il)'^ {t)tpr{t) ■ (b) The free propagators. 

Meanwhile, all vertices in the expansion should be con- 
nected in pairs by the electron propagators. There are 
only two types of the propagators involved: the solid line, 
starting from the vertex qriip^^pi to the vertex q^ftp'^tpr, 
represents the electron propagates in the drain, and the 
dash line, starting from the vertex q^ii^i'i^r to the vertex 

qri4't'^i^ is the propagator in the source. These propa- 
gators connecting two vertices with both the time argu- 
ments s and s' located at the positive (negative) branch 
of the closed time path represent the (anti-) time-ordered 



Green's function G^^^^^(s, s'). The arrow of the prop- 
agator coincides with the propagating direction of the 
tunneling electron. While the propagator connecting two 
vertices, one located at the positive branch and the other 
at the negative branch separately, represents the corre- 
lation Green's function G^^^^^~^{s, s') with respect to the 
arrow pointing to the vertex at the positive (negative) 
branch. The diagrammatic representation of the free 
propagators is summarized in Fig. [T](b). Different from 
the usual closed-time-path contour used in the literature, 
we choose an alternate contour p depicted in Fig. [TJ The 
lower (upper) axis represents the positive (negative) time 
branch of the closed-time-path contour. The contour or- 
dering t >p t' denotes the time t' is former than t along 
the arrow of the contour p, and the corresponding oper- 
ator q(t') is applied earlier than q{t). This prescription 
makes the calculation more convenient in the treatment 
of the time ordering of coupling operators discussed be- 
low. 



B. The real-time diagrammatic expansion 

Now, we can diagrammatically calculate the partial 
trace in the expansion of Eq. PH)) 



Tr, 



(24) 



where H'g_^ ■ ■ ■ H'^^ comes from the perturbation expan- 
sion of the evolution operator Tp{eyip{—i jT^^ dsH^)-} 
with t > T > {si, ■ ■ ■ , Sn)- Conveniently, the vertex with 
the time argument Si=i^...,„ is called the internal vertex 
and that with (i, t) is called the external vertex. The 
time contour ordering Tp in Eq. comprises all the 

permutation of the time series {t > t > (si,-- - ,s„)} 
along the closed time path p, namely, one must sum all 
allowed time contour orderings Tp — Tpj. In terms of 

(jrii'^i'i and q^ii^^i^r, each component of Eq. (|24p can 
be expressed as 



Trs 



Tp[(qri{t)i:tmi{t) 

X (9r';'(T)V',1;(T)V^p(T) 



{+) 

(+) 



X (9r„/„(s„)?/',^ (s„)?/'/„(s„)) Ptot(io)| 



■rnln'i'i'l 



(25) 



where Srir'i'---r„i„;i associated with coupling operators is 
defined as 

(26) 



5 



and the corresponding coefficient £rir'i'---rni„;i is the con- 
tribution of integrating out the electron reservoirs and is 
defined by 



Trf 



(27) 



Note that Sri----i in Eq. ((26)) consists of n + 2 vertices 
which is ordered according to Tp^i. Due to the parti- 
cle number conservation in electron tunneling processes, 
only even orders (n =even) in the perturbation expan- 
sion have contribution. Half of these n + 2 vortices will 
carry with the coupling operator g"*", and the others with 
q. Each Sri--- yiSri--- -i is expressed by a set of topology- 
independent diagrams, in which each diagram are com- 
posed of several allowed closed loops connecting n+2 ver- 
tices. The topology-independence means the order and 
the direction of all propagators lines and the loop assem- 
bly are different. The coefficient £ri----,i can be directly 
calculated only from this set of topology-independent 
diagrams by summing all the topology-independent di- 
agrams with a prefactor (— where I is the 
loop number in the individual topology-independent di- 
agram. The prefactor (— 1)("+2)/2 comes from the factor 



(-0 



(ra-l-2) 



of the n-th order perturbation, and the prefac- 



tor (—1)' is due to the permutation between the fcrmion 
operators in the contraction. An explicit example of cal- 
culating £ri--- -,i'Sri--- ;i for the order of n = 2 can be found 
in Appendix B. 

Accordingly, we introduce a loop operator Cr{ - ■ ■) to 
calculate the total contribution of all time contour order- 
ings in the n-th order perturbation K^^'^\t — r) * Ps[t). 
This loop operator is defined by a loop in particular 
topology-independent diagrams along the real time axis, 

£ I Si S2 S3 ■■■Si--- S2p Si 

fc]^ k-^ k'2 ' ' ' kj ' ' ' kp k/p 



E 

ki---kp k[---k' 



E 



dsi 



ds 



2p 



to 



xTr 



[D[ak, (si)]Gi*;^ (si, S2)Z?Ki (s2)]G(°) {s2, S3) 



■D[ak'(s2p)]G'j^}is2p,si) 



(28) 



where Si=i,...^2p < t are the time arguments of the in- 
ternal vertices, p is an integer with 2p < n, kj and 
k'^ are the energy indexes with (kj ~ Ij , k^ = Tj ) or 

{kj = ^ji^'j — h)^ G'^°''(s,s') is the Keldysh's ma- 
trix shown in Eq. ([TS]), and the 2x2 functional- 
derivative matrix £'[a(s)] is defined by Z?[a(s)] = 

^ \ . 

5 with a{s) being a time-dependent pa- 

Sa*{s) ) 

rameter. The loop involving the external vertices is 

t Si ■ - - T ■ ■ ■ S2p-2 t 

h 



6a{s) 





denoted as £-r 



■ ■ ■kj ■ ■ ■ 



Hp 



This 



definition is the same as Eq. l[28)) but without taking 
the time integral for the time arguments t and r. The 
loop operator is written in such a way that a pair of ver- 
tices labeled by the time arguments s' and s together 
with the propagator with the energy index k represent a 

■ s' s • ■ 



segment of the loop operator Cr 



Each 



internal (external) vertex with the coupling operator qri 
or is described by the operator D[ai] {D[-fi]) and 
£'[ar] (D[7r]), respectively. The connecting propagator 

is the Keldysh's matrix gI^\s, s'). Along the propagat- 
ing direction in each loop, the corresponding loop op- 
erator is written down in the order from right to left. 
The loop is end-point-independent, with which the def- 
inition of the loop operator coincides. Since the hermi- 
tian of the physical quantity, K^^\t — r) * Ps(t) can 
be separated into two parts which are hermitian conju- 
gate to each other. This leads to two sets of the dia- 
grams which are dual each other. The duality of two 
diagrams is defined that the replacements of the vertices 
q {(t') in one diagram by the vertices q^ (q) equals to 
each other. Therefore, we only need to calculate one 
set of the diagrams. In terms of the loop operators, 
only topology-independent diagrams along the real time 
axis should be taken into account. Furthermore, besides 
the prefactor ;^(— l)*-"^^-*/^^', where the factor comes 
from the n-th order perturbation, a weight factor (the 
number of topology-equivalent diagrams) should also be 
added for each topology-independent diagram to calcu- 
late i^(")(t — t) * ps{t) correctly. 

Meanwhile, each order of the kernel expansion 
t)*ps{t) in Eq. (fTT|) contains multi-particle (2, 4, • • • , n+ 
2 particles) correlations. We can re-express the pertur- 
bation expansion of K{t — r) * psiT) in terms of the ir- 
reducible diagrams 



K{t -t)* ps{t) 

K'i\t - r) * ps{r) + K\l\t - r) * ps{t) -f 



(29) 



Explicitly, the n-th order perturbation K[^\t — r) * 
Ps{t) contains n -I- 2 vertices with the time arguments 
{i, r, si, ■ • ■ , s„}. All vertices are denoted by filled cir- 
cles without discriminating q and Also, n internal 
vertices are treated indistinguishably, namely, a diagram 
which exchanges arbitrary two time arguments of inter- 
nal vertices is topology- invariant. The counterclockwise 
and clockwise loops through the same vertices are also 
equivalent except that the vertex orderings along the 
both loop are different. Thus, we define irreducible dia- 
grams as all connected topology-independent diagrams. 
The connected diagram means that each loop of the di- 
agram should intersect with other loops at least once. 
As a result, K[^\t ~ t) * ps{t) comprises all irreducible 
diagrams with the prefactor (— 1)("+2)/2+', Each loop of 
irreducible diagrams is then given by the irreducible loop 
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operator defined as follows: 

i'ir,T {t , Si , ■ ■ ■ T ■ ■ ■ Si • ■ • , S2p-2 , t ) 

t Si ■ ■ ■ T ■ ■ ■ Si ■ ■ ■ S2p-2 t 

h ■■ ■ Tp 



t S2p~2 ■■ - Si - - -T - ■■ Si 



t 



(30) 



The loop operators up to the second order perturbation 
can be easily calculated accordingly, see Appendix C. The 
explicit result of the irreducible diagrams and the corre- 
sponding loop operators is shown in Fig. [2l 



(a) 
(b) 






- ^irj Si , S2 , r, - 4> J- r, Sj , S2 , 

FIG. 2: Irreducible diagrams and the corresponding loop op- 
erators, (a) and (b) The leading and the second order contri- 
butions. 



C. Master equation for the reduced density 
operator 

Thus, writing down the corresponding loop operators 
one by one according to the resulted irreducible diagrams, 
and applying these irreducible loop operators to the fol- 
lowing generating functional 



J(<?,g+;a,7) = y dt'C[7(i')] \ rfr'C[f (r')] 



'{7.Wg.^W(-)-7;W(-)9.^W},(32) 



+ 



and then taking all parameters (a(s),7(s)) to be zero, 
the exphcit expression of K^^^ (t — r) * ps(t) can be ob- 
tained in terms of coupling operators. Here, the func- 
tional derivatives | ^J^^,) • ■ • j^ip; | in irreducible loop op- 
erators are responsible to generate the correct orderings 
of the coupling operators, namely, {Srir'v ■■■r^i^^i)- 

As a result, we obtain the master equation for the re- 
duced density operator expressed in terms of the irre- 
ducible loop operators to all orders in perturbation ex- 
pansions, 

^^Ps{t) = -i{Hs.Ps{t)\~ f dTK{t~T)*ps{T), (34) 



K{t^T)*ps{T) = Qs{t-to) 
x{/C,.(i,T)-e^M. J(q,<z+;a,7)} 



7— 7' — {a}— 



(35) 



where ICir{t,T) consists of all allowed irreducible dia- 
grams (in terms of the irreducible loop operators) 



^„(<,r) =^(-l)"Ar,r (t,Pr{r,Sl,--- ,S2„},t) 
n=0 

00 



n,n'— 



(36) 



with P^{r, si,S2,S3,--- ,Sn} = {t, si,S2,S3,---} -I- 
{si, T, S2, S3, ■ ■ ■ } + ■ ■ ■ containing n + 1 permutations, 
and the operator VV(t) is defined as 



w(<) = 



(-1) 



n+l 



[si,---S2„,si], (37) 



n=l 



a log form in the Taylor expansion, and the loop operator 
l^in,T \ • •] only involving the internal vertices is given by 





c 


r[si ,• 


■ ■ S2p, Si ] 








S2p, Si ) . 






r,T ( -^l 




= tr 


(si 


S2 • • • S2p Si \ 


(31) 


\ ki 


K ) 




= Cr 




S2p ■ ■■ S2 Si \ 






ki J 



(38) 



In fact, the operator e^^*-* generates all the loop opera- 
tors for the reduced density operator ps{t): 



Ps{t) = Qs(i-io)|e^(*)j™((7,9+;a) ^ ^ }. 

L {a}=0J 



(39) 
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The leading order contribution (n = 0) to the master 
equation is obtained as foUows 



d_ 
dt 



ps{t) - -I [Hs, ps{t)] -Qs{t- to) 



f dTKt\t-T,Ps{T)l (40) 

J to 



(0) 



: I dT'C[f{T')]Qs\r-to)ps{r)} 



7— 7'— 

(41) 

Also, the time variated reservoir fluctuation due to the 
interaction with the qubit has been taken into account. 
The internal vertices of ICir{t,T) and e^*^'^^ in Eq. ([55]) 
are mixed together for higher order contributions. The 
reduced density operator can not be extracted unless the 
following approximation is utilized, 



7— 7' — {a}— 



7— 7' — {a}— 



(42) 



Jsa('?,9+;«,7) = J dt'Cmt')] J dr'C[f(r')] 

X Te^' (''^'Jw ('')1" 

X {QsHr~to)psir)] 

(43) 

Without resorting the traditional diagrammatic tech- 
nique in the Laplace space, the real-time diagrammatic 
technique has been developed to derive the master equa- 
tion. The charge qubit dynamics in the non-Markovian 
regime can be studied based on Eqs. (|34j|4ip . 



IV. QUBIT DECOHERENCE 

To explore the qubit decoherence induced by the QPC 
measurement, the charge qubit as a single electron in a 
double quantum dots[8| is considered. The Hamiltonian 
of the system in Eq. ([2]) can be explicitly written as 

Hs = Elc+cl + Ebc+cr + no{c+CR + c+cl), (44) 

where cJ(cl) and cJ(c_r) are the creation (annihilation) 
operators of the electron sited in the two dots labeled 
by L and R with the single-electron constraint cj^CL + 



CrCr 



1, El,r are the corresponding energies, and VLq is 
the electron hopping amplitude between the double dots. 
The interaction between the system and the QPC due 
to the measurement is characterized by the interaction 
couphngQ 



(45) 



and CI and — il' in Eq. (|45|) are the electron hoping 
amplitude of the QPC without and with the measure- 
ment of the single electron in the double dots. Eq. (|45|l 
describes a variation in the barrier of the QPC when the 
single electron occupies on the right dot. 

We shall consider the leading order contribution to the 
master equation. According to Eqs. (|40l4ip . we obtain 

^os{t)^-i[Hs,ps{t)]- f dT\Ro,[[k{t-T) 



dt 



X R{t-T),e''"'^'-^^psiT)e 



iHs{t-T) 



]] , (46) 



where the double bracket [[A,B]] ^ AB - {AB)+ , the 
operator R[t) is given by 



R{t) = cos0(|e)(e|-|5)(ff|) 
s^^* \g) {e ■ 



— sm ( 



e-^^* |e) {g\) 



(47) 



with 7 — ^4r2§ + {El — ErY the energy difference be- 
tween the ground state \g) and the excited state |e) of 
the qubit, 9 = cos-^[{El - Er)/-^], and i?o = R{t = 0). 
The reservoir correlation function k{s) in Eqs. (|46p which 
characterizes the QPC structure associated with the tem- 
perature effect and the external bias is expressed as 



dede' 
1^' 



,i(e-e')s 



fc(e,e'), 



(48) 



fRm-fL{e'))J(e,e') (49) 
+fL{e){l- fR{e'))J{e',e), 



where fc(e, e') is the electron-tunneling spectrum for the 
QPC, fLM{() = 1/(1 + cxp/?(e - p.l,r)) are respectively 
the Fermi-Dirac distribution functions for the source and 
the drain, and the spectral density J(e', e) for the QPC 
structure is given by 



J{e,e')^i:^gL{e')gR{t)\^'{e,e')\ 



(50) 



with gL.Fti^) being the density of states of the source 
and the drain. Also, we assume here the electron energy 
levels {ei,er} are continuous. Eq. ((46l) contains the non- 
Markovian processes of the qubit dynamics up to the 
leading order. It can be checked Eq. reduces to the 
result in Refs. [l^, [l3| in the Markovian limit. 
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FIG. 3: (color online). The qubit dynamics and the corresponding reservoir time correlation function in the measurement 
condition with an extreme small ui (w ^ 7) and different a. The symmetric coupled quantum dots (7 — 2flo) is simulated. 
The qubit is set initially in the \L) state. The following measurement parameters are used with Qq being a rescaling factor: 
/3 = l/f^o, Vd = 205^0 and p = Qq- (a) and (b) The qubit dephasing. The time evolution of the real and imaginary parts of the 
qubit density matrix element peg are plotted, respectively, (c) The qubit relaxation. The time evolution of the qubit density 
matrix element pgg for the qubit in the ground state is plotted, (d) The amplitude of the reservoir time correlation function is 
plotted in units Qq. (e) The time derivative of the phase in the reservoir time correlation function is plotted in units Qo- The 
integration in Eq. (|56p is approximately calculated by fitting the Lorentzian spectral density with Gaussian function to remove 
the energy cut-off. 



To make the non-Markovian feature apparent, we shall 
concentrate on the charge qubit with symmetric coupled 
dots El = Eji characterized by 7 = 2ilo {0 = 7i'/2). 
The equation of motion for the reduced density matrix 
becomes. 



dt 



'i"fPeg{t)~2 / dT\k(t — t)\ COS (l){t — t) 

J to 

(51) 



dt 



2 / dT\k{t-T)\ 
Jto 

I COs[(t - r)7 - (j){t - T)]pee{T) 

cosiit - t)j + cj^it - T)]pggiT)Y (52) 



where the matrix elements are defines as pij {t) — 
{i\ps{t)\j) with being the ground state or the excited 
state of the qubit, and <j>{t) is the phase of the reservoir 
correlation function k{t) = |fc(t)|e~*'^(*^ 

The qubit decoherence can be studied by the analy- 
sis of the spectral density J{er,ei). In the Literature, 
the density of states in the QPC reservoirs and the hop- 
ping amplitude across the QPC barrier are assumed to be 



energy-level independent, namely, the wide-band approx- 
imation for the QPC structure. [1, E [H [11 El Q 
The spectral density is then given by 

J{er,ei)^7r^gL9R\n'\\ (53) 

This corresponds to the Markovian limit, in which the 
qubit dynamics is strongly decoherent [Til. [H [H, [l^ . 
The Markovian dynamics arises from tunneling electrons 
in the QPC with a shortest time correlation k{t — t) oc 

S{t-T). 

However, as indicated in Eqs. (|48l49p . the QPC struc- 
ture determines the correlation time scale of the tunnel- 
ing electron fluctuation in the QPC. The non-Markovian 
processes of the qubit dynamics emerges only when the 
QPC structure is casted with a finite correlation time 
scale. The effect of the QPC structure can be character- 
ized by an energy-level dependence of the spectral den- 
sity. In the literature, the spectral density Jrd{^) for a 
reservoir coupling to a dot (or a molecular wire) is pa- 
rameterized by Lorentzian spectrums.[l7j 



PkCtk 



(54) 



where pk, Wk and are fitting parameters. The 
Lorentzian spectrum has also been applied to study 



the quantum measurement of this system (using a con- 
stant hopping amplitude with a Lorentzian density of 
state). [2^ For our system which involves electrons tun- 
neling a barrier between two reservoirs, the spectral den- 
sity can be approximately treated as 



pa 



2 {\er-ei\-wy + ay4: 



(55) 



In Eq. (|55p . the parameter p specifies the magnitude of 
the spectral density, w characterizes the variation of the 
barrier potential due to the interaction with the qubit 
electron, and the width a indicates a modulation of the 
decay rate for the qubit decoherence. Changing the 
variables of the integration in Eq. from (e,e') to 

(e' — ^^,£ — ^Y^) and integrating out e', we then ob- 
tain the following reservoir correlation function. 



k{s) = 



de 
— f 
27r 



pa 



4(|e| - w)2 -f a2 



:[g{e-Vd)+gie + Vd)}, (56) 



where the function g{x) = ■ 

A close connection between the qubit decoherence and 
the tunneling-electron fluctuation is revealed in Eq. ([56| . 
According to the spectral density in Eq. ([55]) . the ran- 
domness of tunneling electrons across the QPC barrier is 
expounded first. Obviously, taking the limit a <C je; — e^l 



to Eq. ((55l) . J{er,ei) reduces to tt S{\ei ~ e^l — w), only 
two channels e; — = involved. On the other hand, 
a \ei — er \ leads to a channel-mixture regime, where all 
transitions e; ^ that electrons tunneling between the 
source and the drain are allowed with the weight deter- 
mined by J{er,ei). The randomness of electron tunnel- 
ing processes in the channel-mixture regime comes from 
electron scattering which are determined by the band 
structure associated with the geometry of the two metal 
gates in the QPC and modulated by the interaction with 
the qubit electron. The parameter a in Eq. (j55p char- 
acterizes the deviation describing how transfer energies 
je; — e^l in all electron tunneling processes are close to w 
in the statistics of ensemble average. The qubit dynamics 
with different a and an extreme small w{<^ 7) are simu- 
lated in Fig. [3] with (a) and (b) for the qubit dephasing 
and (c) the qubit relaxation. The result shows that the 
qubit decoherence is suppressed as tunneling electrons 
with a smaller deviation. Associated with the qubit de- 
cay rate, a can be used to judge the decoherent behavior 
of the qubit state. 

Furthermore, the reservoir time correlation function 
in Eq. (j56|) describes the time correlation of timneling- 
electron fluctuation in the variation due to the measure- 
ment. In Fig. 21 the reservoir time correlation function 
is plotted in the condition that the variation of the QPC 
barrier potential has the same energy scale as the qubit 
system. It can be found in Fig. 2] (a) that the amplitude 
|fc(s)| is a Lorentzian profile with periodic deep peaks at 
t = tt/2w, 3tt/2w, The amplitude \k{t)\ describes 




20 



is 



' 



t (in units : 1/Q ) 

FIG. 4: (color online). The reservoir time correlation function 
plotted in the measurement condition that the variation of 
the QPC barrier potential has the same energy scale as the 
qubit system {w = 7) and with different a. The symmetric 
coupled quantum dots (7 = 200 ) is considered. The following 
measurement parameters are used with Qq being a rescaling 
factor: (3 = l/Oo, Vd = 20r2o and p = Oo, respectively, (a) 
The amplitude of the reservoir time correlation function is 
plotted in units Qq. (b) The time derivative of the phase of 
the reservoir time correlation function is plotted in units flo. 



the correlation time scale between tunneling electrons 
which are induced by the measurement in the time in- 
terval t. With fixing the amplitude, the wider the half 
width of the Lorentzian profile of \k{t)\, the longer the 
time correlation of tunneling-electron fluctuation. Fig. [4] 
(a) then indicates that the time correlation increases with 
a decreasing. In other words, a shorter time correlation 
between electron tunneling processes leads to a more ran- 
dom electron-tunneling spectrum (with wide profile), see 
Eqs. (|48|49p . In the channel- mixture regime, the time 
correlation of tunneling-electron fiuctuation is smeared, 
and the reservoir memory effect on the qubit dynamics 
is suppressed. The qubit dynamics with different a are 
simulated in Fig. [5l The result shows that in the channel- 
mixture regime (a — 2O7), the qubit undergoes a severe 
decoherence, which corresponds to the qubit dynamics in 
the Markovian limit. On the other hand, instructively, 
if the QPC structure can be recasted such that the mea- 
surement approaches to the case almost without random 
electron scatterings {a = O.57), the non-Markovian pro- 
cesses of the qubit dynamics emerges. The qubit simply 
performs a periodic oscillation with small fluctuations. 
The qubit decoherence is suppressed. The mechanism of 
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the qubit decoherence can be understood by the analysis 
of the time derivative of the phase (f>(t) in the reservoir 
time correlation function, which is plotted in Fig. [4] (b). 
According to Eq. ([5^ . d(j){t)/dt characterizes the aver- 
age transfer energy that the qubit electron obtains from 
all tunneling processes in the QPC. Fig. |4] (b) indicates 
that in the channel- mixture regime (a = 2O7), the neg- 
ative energy-transfer rate before t = 0.25/r2o leads to 
the qubit electron energy being exhausted. The contri- 
bution that the qubit electron absorbs energy from the 
QPC is very small and ignorable after t = 0.25/r2o> be- 
cause \k{t)\ decay to zero, see Fig. |3] (a). This effect 
occurs repeatedly due to the non-Markovian dynamical 



structure. The qubit finally relaxes to a mixed state, 
as shown in Fig. [5l However, for the case a = O.57, 
the qubit electron absorbs (emits) energy from (to) the 
QPC before (after) t = 7r/2w, 37r/2w, • • • in each period, 
as shown by the peaks with positive (negative) energy- 
transfer rates in Fig. |4] (b). The qubit is driven by the 
QPC periodically. Because of a longer time correlation 
of the tunneling-electron fluctuation [see the wide pro- 
file of \k{t)\ in Fig. [4] (a)], this periodic driving becomes 
significant. The qubit state is oscillating almost without 
decoherence, if the effect of the random electron scatter- 
ings can be ignored. 




Markovian regiine: 

a= 20y 

non-Markovian regiine: 

a= Y 

a= 0.5 y 



t (in units : 1/D ) 30 




FIG. 5: (color online). The qubit dynamics measured in the condition = 7 and with different a. The symmetric coupled 
quantum dots (7 = 2Qq) is simulated. The qubit is set initially in the |L) state. The following measurement parameters are 
used with Qq being a rescaling factor: (3 = l/Jlo, Vd ~ 20f2o and p — flo, respectively, (a) and (b) for the qubit dephasing, and 
(c) the qubit relaxation. 



Accordingly, the qubit decoherence also depends on the 
number of the peaks, which is determined by the varia- 
tion of the barrier potential. By comparing Fig. [HKd) and 
(e) with Fig. 21 it can be found that the number of the 
periodic peaks inside the half-width of the reservoir time 
correlation function indeed decreases with the parameter 
w decreasing. The more the sharp peaks with positive 
energy-transfer rates that the reservoir time correlation 
function contains, the larger amplitude of the oscillation 
the qubit performs in the asymptotic regime, specially, 
the population Pgg(t) and the coherence Iie[peg(t)]. A 
large amount of energy transferred between the qubit 
electron and the tunneling electrons in the case with 
large w leads to a severer time variation (oscillation) on 
the qubit-electron population and the induced coherence. 
However, the w-dependence of the qubit dynamics is not 
so sensitive for the large a. The qubit decoherence under 



the QPC measurement with different w almost coincide, 
as indicated by the case of a = 2O7 in Fig. [3] and Fig. 
[5l To drive the qubit state efficiently, tunneling electrons 
must be in the state with a longer time correlation, no 
matter how large the variation of the barrier is. Also, 
Fig- El (b) demonstrates that a sinusoidal quantum oscil- 
lation of the qubit coherence is still accessible under the 
quantum measurement, if the QPC structure is casted 
with a small w of the spectral density. 

Finally, the numerical result also shows that the larger 
the bias voltage is applied, the shorter the time corre- 
lation of the tunneling-electron fluctuation becomes. A 
large amount of electrons tunneling across the QPC bar- 
rier forced by the large bias voltage leads to the time 
correlation of the tunneling-electron fluctuation being 
smeared with a sharp proflle |fc(i)l- Therefore, the qubit 
dynamics becomes typically Markovian, and the qubit 
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state is severely decoherent. The same decoherence be- 
havior also occurs at a higher temperature. As a re- 
sult, the constant spectral density used in literatures 
essentially describes electron tunneling processes with 
the most random channel mixture, and the qubit sim- 
ply undergoes decoherent processes. More realistic non- 
Markovian dynamics emerges only when the effect of the 
explicit spectral density is taken into account. A particu- 
lar design to cast the QPC structure such that tunneling 
electrons are in the state with longer time correlation 
could suppress the qubit decoherence during the QPC 
measurement. 
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APPENDIX A: DERIVATION OF EQUATION OF 
MOTION IN EQUATION ([7]) 



SUMMARY 



We have developed a non-equilibrium theory for the 
charge qubit dynamics accompanied with the QPC mea- 
surement. The effect of non-equilibrium fluctuation of 
the QPC reservoirs to the qubit dynamics has been 
treated by using the real-time diagrammatic technique 
developed based on Schwinger-Keldysh's approach. We 
introduce the loop operator to exactly derive the master 
equation for the reduced density operator, which is ex- 
pressed in terms of the irreducible diagrams. The qubit 
decoherence is studied according to the resulted master 
equation up to the leading order. The non-Mar kovian 
processes in the qubit dynamics has been taken into ac- 
count in this framework. We find that the qubit deco- 
herence sensitively depends on the spectral density of the 
QPC structure. The constant spectral density in fact 
describes the electron tunneling processes with the most 
random channel mixture, which causes the severest qubit 
decoherence, and corresponds to the qubit dynamics in 
the Markovian limit. However, if the QPC structure can 
be controlled such that the spectral density with a nar- 
rower width, the electron tunneling processes reduce to 
a less channel mixture. The non-Markovian processes 
of the qubit dynamics emerges. The qubit simply per- 
forms a periodic oscillation with less decoherence. A 
longer time correlation of the tunneling-electron fluctua- 
tion results in the qubit state being periodically driven. 
In the channel- mixture regime, the time correlation of 
the tunneling-electron fluctuation is smeared by random 
electron scatterings. The qubit fails to be driven effec- 
tively by the QPC due to short time correlations of the 
tunneling-electron fluctuation, and the qubit state is sim- 
ply decoherent. As the measurement operation is de- 
signed with the minimum deviation of electron tunneling 
processes, the qubit state could be measured with the 
least decoherence effect. The further work taking into 
account higher order contributions to the qubit decoher- 
ence and the noise spectrum of the QPC output signal is 
in progress. 



The equation of motion for the reduced density oper- 
ator listed in Eq. ([7]) is derived in this appendix. In 
the interaction picture, the formal solution of the Liou- 
ville equation ■§iPtot = —i[H,ptot\ for the total density 
operator is given by 



Ptotit) = Ptotita) - i / dT[H'^,ptot(T)], 



(Al) 



iteratively, which leads to Dyson's perturbation expan- 
sion for the total density operator. However, to system- 
atically build up the perturbation scheme for the kernel 
expansion in terms of the reservoir contour-order Green's 
functions (two-point correlation functions), alternatively, 
the second iteration of Eq. (jAl[) 



Ptot{t) = Ptot{to) - i / dT[H!^,ptot{to)] 

Jto 

dT r dT'[H',,[H',,,ptot{T')]], (A2) 

Jto 

which contains Hamiltonian operator orderings of H'^ and 
H'^i in different time should be used. Obviously, as indi- 
cated by Eq. (|A2p . the total density operator obeys the 
following equation of motion 



d_ 

dt 



Ptotit) 



H'nPtotih) 



dr 



to 



KrPtotir) 



(A3) 

Then, taking the partial trace to integrate out the degrees 
of freedom of the QPC reservoirs, Eq. ([7]) is resulted. 



APPENDIX B: EXAMPLE OF CALCULATING 
THE COEFFICIENT Sh-.-a 

In this appendix, we present an example of the contri- 
bution of the perturbation expansion Eq. (j25p for n — 2 
with the particular contour ordering {t >p r >p to >p 
{sl^ S2)} to illustrate how to diagrammatically calculate 
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the coefScient Srir'i'rihr2i2 Eq. ([27]), 

Srlr'l'rilir2l2 —P ^0 i^l: ^2)^ 

xSrlr'l'rihrihi^t >p T >p tf) >p (si,S2)^ 

= '?HW'7r-'i'(T)ps(io)r|g+;^ (51)9^2(2 (S2)} 



xf { (si)) (V3+ (S2)4 (S2)) }] . (Bl) 

The corresponding diagrams are shown in Fig. [51 It in- 





FIG. 6: Four topology-independent diagrams of the example 
Eq. ([Bit . 



eludes four topology- independent diagrams, two single- 
loop diagrams with the prefactor (— l)(2+2)/2+i = 
— 1 and two double-loop diagrams with the prefactor 
(_l)(2+2)/2+2 _ ^ rpj^g explicit expression of the co- 
efficient £rir'i'rihr2i2 Eq. (jBl[) that the partial trace 
is carried out can be directly written down from four 
topology-independent diagrams 

£rlr'l'rilir2l2 *0 >p (si, S2)^ 



xG^°);Js2,i)G^^;,(i,s2) 

-<5;i/2<^r2r<5/i''5rir'G^^^^ (s2, Si)G<_V2 (i, S2) 



(0) 



(0) 



Shl'Srr'6rir2Sl2lG''^\> (t, Si)G^]., {t, t) 

(0) +\r^(o) 



xG^",^,^(s2,i)G^^; (S1,S2). 



(B2) 



All coefficients Srir'i'--- for higher order contributions 
with arbitrary contour ordering can be calculated dia- 
grammatically in the similar way. 



APPENDIX C: THE LOOP OPERATORS FOR 
LEADING AND SECOND ORDER 
PERTURBATIONS 

As an example, we illustrate how to write down the 
loop operators for the leading order kernel K'^^^ {t ~ t) * 
Ps{t) and the second order kernel K'^^^t — r) * ps{t) 
diagrammatically. For the leading order kernel, only two 
topology-independent diagrams involves, and one is dual 
to the other. One of the both is shown in Fig. [7] (a). 



(a) 



IrT 



0! 



(b) 




2 %w\u/ '\«?/' 

• " • o o 

2? T(sf 




, 2 %iJ'~X.u H'-^^'V 

^- ^' /, 





FIG. 7: Topology-independent diagrams of K{t — r) * ps(j) 
up to the second order perturbation, (a) and (b) The leading 
and the second order contributions 

For the second order kernel with all possible time con- 
tour orderings, it contains totally 24 diagrams, half of 
them are dual to each other. We only need to consider 
one set of them (12 diagrams). Furthermore, one can 
find that only six topology-independent diagrams exist, 
see Fig. [7] (b). Each one corresponds to two topology- 
equivalent diagrams, i.e. the weight factor is 2 for each 
topology-independent diagram, while the prefactor is ^ 
for the single-loop diagrams and if for the double-loop 
diagrams. 

The loop operators can be directly written down ac- 
cording to the topology-independent diagrams in Fig. [T] 
For the leading order kernel, the corresponding loop op- 
erator is given by ^ ^ ^ y ^ ^ . We thus obtain 
the following expression 

K^^\t - t) ^ ps{r) 
xj(g,g+;a,7) ^ ^ 

7— 7' — {a}— 



T t 

r' V 



dualityj. (CI) 



Similarly, the explicit expression of the second order 
kernel can also be obtained according to the topology- 
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independent diagrams in Fig. [7| (b), 
K^^\t - t) * Ps{t) 

+duality| , 
21 \ n h \ r2 h 



7— 7' — {a}— 

(C2) 



-I- 



t S\ T S2 t 

ri h r2 h 



32 Si S2\fi t T t 

n h r^'[ r' I' 



t Si S2 T t 

n h r' I' 



T S2 
I' 



t T S2 Si t 

r' V ri h 



Note that all contour-ordering {Tp.} of the leading and 
the second order perturbations have been taken into ac- 
count in Eqs. (|C1IC3|) combining with their duality, re- 
spectively. 

In addition, we can express each order of the kernel 
in terms of the loop operators associated with the irre- 
ducible diagrams, for instance, 

7— 7' — {a}— 



{t, T, t) 

-Cir^r {t, Sl, r, S2, t) - Cir^r (t, T, Si, S2, t) | 



x>7(9,9^;<3,7) 



7— 7' — {a}— 



(C4) 



(C3) 



Therefore, we can resum all perturbation orders of K{t- 
t) * psiT) according to the irreducible diagrams. 
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